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Abstract 

Background: Quantitative measurements in the myocardium may be used to detect both focal and diffuse disease 
processes that result in an elevation of T1 and/or extracellular volume (ECV) fraction. Detection of abnormal 
myocardial tissue by these methods is affected by both the accuracy and precision. The sensitivity for detecting 
abnormal elevation of T1 and ECV is limited by the precision of T1 estimates which is a function of the number and 
timing of measurements along the T1 -inversion recovery curve, the signal-to-noise ratio (SNR), the tissue T1, and 
the method of fitting. 

Methods: The standard deviation (SD) of T1 and ECV estimates are formulated and SD maps are calculated on a 
pixel-wise basis using the Modified Look-Locker Inversion recovery (MOLLI) method. SD estimates are validated by 
numerical simulation using Monte-Carlo analysis and with phantoms using repeated trials. SD estimates are 
provided for pre- and post-contrast optimized protocols for a range of T1s and SNRs. In-vivo examples are provide 
for normal, myocarditis, and HCM in human subjects. The formulation of SD maps was extended to R1 and ECV. 

Results: The measured myocardial SNR ranged from 23 to 43 across the heart using the specific T1 -mapping 
protocol in this study. In this range of SNRs, the estimated SD for T1 was approximately 20-45 ms for pre-contrast 
myocardial T1 around 1000 ms, and was approximately 10-20 ms for post contrast T1 around 400 ms. The proposed 
estimate of SD was an unbiased estimate of the standard deviation of T1 validated by numerical simulation and 
had > 99% correlation with phantom measurements. The measured SD maps exhibited variation across the heart 
due to drop off in surface coil sensitivity as expected for the variation in SNR. Focal elevation in T1 and ECV was 
shown to have statistical significance on a pixel-wise basis for in-vivo examples. 

Conclusions: Pixel-wise estimates of T1 mapping errors have been formulated and validated, and the 
formulation has been extended to ECV. The ability to quantify the measurement error has potential to determine 
the statistical significance of subtle abnormalities that arise due to diffuse disease processes involving fibrosis 
and/or edema and is useful both as a confidence metric for overall quality, and in optimization and comparison 
of imaging protocols. 
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Background 

Quantitative methods such as Tl -mapping and extra- 
cellular volume (ECV) mapping appear promising to 
complement LGE imaging in cases of more homoge- 
neously diffuse disease which affect the myocardial 
extracellular space [1-17]. Quantitative measurements 
in the myocardium may be used to detect both focal 
and diffuse disease processes that result in an eleva- 
tion of Tl and/or ECV. 

Direct measurement of extracellular volume (ECV) was 
initially developed for quantifying the myocardial extra- 
cellular fractional distribution volume [1] and has been 
proposed as a means for detection and quantification of 
diffuse myocardial fibrosis [2-6,15-17]. This approach is 
based on the change in Tl following administration of an 
extracellular contrast agent and circumvents the limitation 
of LGE in detecting a global change in Tl, which typically 
uses a single post-contrast Tl measurement. The myocar- 
dial ECV is measured as the percent of tissue comprised 
of extracellular space, which is a physiologically intuitive 
unit of measurement. ECV has been shown to correlate 
with collagen volume fraction [3,5]. The topic of ECV 
mapping as well as native Tl -mapping is of current inter- 
est as a diagnostic tool for a wide range of cardiomyopa- 
thies including dilated cardiomyopathy [2], diffuse fibrosis 
associated with myocardial dysfunction in congenital heart 
disease [4], aortic stenosis and HCM [5], cardiac amyloid- 
osis [15], Anderson Fabry [17], myocardial infarction 
[8,11], as well as aging processes [8]. 

The sensitivity for detecting abnormal elevation of Tl 
and ECV is limited by the precision of Tl estimates 
which is a function of the number and timing of mea- 
surements along the Tl-inversion recovery curve, the 
signal-to-noise ratio (SNR), the tissue Tl, and the 
method of fitting. Detection of abnormal myocardial 
tissue by these methods is affected by both the accuracy 
and precision. In this work, we only consider the ran- 
dom component due to noise which limits precision and 
not bias errors that affects accuracy. Although absolute 
accuracy of in-vivo measurements is an important and 
open subject, these methods have been shown to be 
highly reproducible [18] in practice despite bias errors. 

Quantifying the statistical fluctuation of measured Tl 
could improve confidence in assessing the significance 
of results. We propose to produce a map calibrated in 
Tl units that represents the quality of the Tl estimate, 
by transforming the standard deviation (SD) of residual 
fitting error into the SD of the estimated parameters. 
The estimate of the Tl parameter error based on the fit 
residuals is derived analytically, and validated on phan- 
tom measurements. Robust methods for both fitting and 
estimating the SD of the underlying measurements are 
described. In-vivo examples are shown to demonstrate 
the potential utility. 



Methods 

Theory 

The sensitivity for detecting abnormal elevation of Tl is 
limited by the precision of Tl estimates, which in the 
case of inversion recovery methods is a function of the 
number and timing of measurements along the Tl- 
recovery curve, the signal-to-noise ratio, tissue Tl, 
method of fitting, and the accuracy of the model. Physio- 
logic fluctuation is not considered. It is proposed to pro- 
duce a standard deviation (SD) map calibrated in Tl 
units that represented the standard deviation of the Tl 
estimate, by transforming the SD of the residual fitting 
error into the SD of the estimated parameters. Estimate 
of the Tl parameter error based on fit residuals is 
derived analytically, and validated by both Monte-Carlo 
numerical simulation and phantom measurements using 
repeated trials. 

Inversion recovery is widely used for Tl -mapping 
using Look-Locker methods. In applications such as 
cardiac MR, a modified Look-Locker (MOLLI) method 
[9,10] uses inversion recovery with multiple single shot 
images at different inversion times. Pixel-wise paramet- 
ric mapping is accomplished by performing a curve fit to 
the multiple inversion time measurements. The original 
MOLLI paper [9] assumes a 3-parameter model of the 
form S(ti) = A - B exp(-ti/Tl*), where Tl* < Tl repre- 
sents the apparent Tl which is shortened by the influ- 
ence of imaging RF pulses. The desired Tl is then 
calculated at each pixel using Tl = T1*-(B/A-1), referred 
to as the Look-Locker correction, originally derived from 
considering a continuous fast low angle shot (FLASH) 
gradient echo readout [19]. The Look-Locker correction 
is used in MOLLI despite the fact that imaging uses 
non-continuous balanced steady-state free precession 
(bSSFP), which violates the assumption used in the 
formulation. This assumption is a key source of bias 
error not treated in this paper and is sensitive to vari- 
ables such as Tl and T2 as described in the literature 
[9,20,21]. In this work, we only consider the random 
component due to noise which limits precision and not 
bias errors that affect accuracy (Figure 1). 

A phase sensitive (PSIR) reconstruction was used [22] 
to restore the sign and thereby avoid performing a 
magnitude fit, or performing multiple fits to estimate 
the zero crossing [9]. In PSIR reconstruction, the real 
component is used which results in normally distributed 
noise [23]. The 3-parameter model may be written as: 

rft)=A-B.exp(-^) =A-B.exp(-^1) (1) 

and the measurements y(t) may be fit for the unknown 
parameters (A, B, Tl) at each pixel using a downhill sim- 
plex minimization approach (Nelder-Mead method) [24]. 
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Figure 1 Illustration of accuracy versus precision. 



The desired covariance matrix C of the estimated 
parameters (A, B, Tl) may be approximated as [25]: 
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where the matrix D 1/2 is a first order approximation to 
the Hessian matrix D. D is comprised of second order 
derivates of the chi-square cost to the parameters. D 1/2 
is computed as: 
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(3) 

after having dropped the second order terms in D. D 1/2 
is comprised of the partial derivatives of the signal 
Eq (1) relative to the estimated parameters, derived ana- 
lytically for each inversion time t { as: 
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The noise of the PSIR reconstructed measurements is 
normally distributed with standard deviation c>i for each 
measurement which can be treated as equal across 



measurements (o*i = a) and are estimated at each pixel 
from the fit residuals. Although the measurements are 
very well approximated as normal, the parameter esti- 
mates themselves may not be normally distributed due 
to the non-linearity of the model. However, in our ex- 
perience the parameter estimates (A, B, and Tl) were 
well approximated as normal in the range of SNR 
encountered. 

The above formulation may also be extended to ECV 
maps, where ECV is calculated as: 



ECV = (1-HCT) 



(R Ipost R lpre ) myocardium 
(Rl^osf-RV-e) hlood 



(5) 



In this case, it is assumed that the random noise com- 
ponent in the ECV measurement is dominated by the 
myocardium since the Rl in the blood is measured in a 
large ROI, which reduces the noise. The calculation of 
SD is reformulated for Rl = 1/T1 by substituting 5y/5Rl 
for 5y/5Tl in Eq (3) above. 
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The SD for ECV is calculated as: 
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which is dominated by the post-contrast variation, 
i.e. a Rlpost ., due to the increased value of Rl following 
administration of contrast. 

Robust estimation 

Robust estimation was used for both fitting for Tl as 
well as for estimation of the standard deviation, a, of the 
measurements. Robust fitting [26,27] uses iterative re- 
weighting to improve the fit in the presence of outliers. 
At each iteration, the weighting of outliers is reduced 
based on the value of their residuals. For instance, out- 
liers may result from errors such as residual uncorrected 
motion. Additionally, a robust estimation of the under- 
lying standard deviation (a) from the fit residuals is used 
based on the median absolute deviation (MAD) ap- 
proach [26,28]. The MAD estimate o is calculated from 
the residuals ^ = (fit - meas) as: 



a = median(abs(ri))/0.6745 



(8) 



where ri are the residuals 8i after discarding the (p-1) 
values with lowest magnitude, p = 3 is the number of pa- 
rameters being fit, and the scale factor 0.6745 is used for 
noise which is normally distributed. The scale factor 
is calculated as the median of the absolute value 
of normally distributed noise with standard deviation 
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equal to 1. Discarding the lowest residuals is necessary 
to avoid a bias error due to over fitting. 

The estimated SD for the corresponding Tl-map value 
is calculated at each pixel producing a SD map, follow- 
ing the flow chart in Figure 2. 

Numerical validation 

A Monte-Carlo simulation using N = 65,536 trials was 
used to compute the standard deviation in Tl (o T1 ) as a 
function of SNR, and Tl for a specific MOLLI protocol 
(5-3 sampling, TImin = 105 ms, Tl shift = 80 ms) and 
T2 = 45 ms, and was compared with the estimate of 
standard deviation based on the proposed approach 
using the fit residuals. Each trial used independent noise 
with a normal distribution. In order to test the robust- 
ness to outliers, a separate test was conducted which 
compared the robust iteratively re-weighted and stand- 
ard methods in the presence of outliers. For this com- 
parison, a single outlier with a randomly positioned 
inversion time (i.e., 1 of the 8 measurements) had a 
noise standard deviation corresponding to 6 times the 
noise standard deviation of the other measurements. 

In addition to measuring the accuracy of the esti- 
mate of SD, i.e. whether the estimate is unbiased, the 
variability or standard deviation of the standard devi- 
ation estimate was measured. 



Phantom validation 

Experimental validation was performed using phantoms 
by repeated measurements (n = 200 repetitions) of a set 
of Tl -measurements comparing estimated and calcu- 
lated standard deviations. A "measured" SD was calcu- 
lated for each pixel, and compared with the mean of the 
estimated SD map values for the corresponding pixel. 
Phantom validation used a set of CuS04 doped agar gel 
phantoms with varying concentrations with Tl and T2 
in the expected range for myocardium, both native and 
with Gd contrast. Phantoms had Tl in the range 
250-1600 ms and T2 in the range 40-75 ms. 

Imaging 

Imaging was performed on 32 channel 1.5 T Siemens 
Avanto and Espree scanners (Siemens Medical Solutions, 
Erlangen, Germany), equipped with 45 mT/m and 
200 T/m/s, and 33 mT/m and 170 T/m/s gradient sys- 
tems, respectively. The MOLLI imaging protocol used in 
this study acquired data with 2 protocols (sampling 
strategies). The first protocol was optimized for native 
(pre-contrast) myocardial Tl values and the second 
protocol was optimized for shorter Tl corresponding to 
Gd contrast. The MOLLI imaging protocol used for 
native Tl acquired data at 8 inversion times over an 11 
heart beat breath-hold at end-expiration with 2 inversions. 
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Figure 2 Pixel-wise calculation of SD values based on fit residuals. 
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For the sake of reduced heart rate variability the 
protocol used a 5(3)3 scheme, acquiring 5 images after 
the first inversion, a 3 second pause for recovery, and 
3 images acquired after the second inversion. The 
protocol used with Gd contrast had a 4(1)3(1)2 sam- 
pling scheme which acquired 9 images in 11 heart- 
beats with 3 inversions which was found to improve 
the Tl estimate. For the contrast protocol with shor- 
ter Tl, a shorter recovery period (Is) is used between 
inversions, which permits additional samples of the in- 
version recovery at short inversion times in the same 
overall breath-hold duration. 

Typical imaging parameters for both sampling schemes 
were: non-selective adiabatic inversion pulse, steady state 
free precession single shot read out with 35° excitation flip 
angle, typical field of view 360 x 270 mm 2 , slice thickness 
6 mm, minimum inversion time 110 ms, inversion time 
increment 80 ms, 7/8 partial Fourier plus parallel imaging 
factor 2. A 2.56 ms tan/tanh adiabatic inversion pulse 
was used to improve the inversion efficiency at the 
low values of myocardial T2. For heart rates below 
90 bpm we typically used: matrix 256 x 144, voxel size 
1.4 x 1.9 x 6.0 mm 3 , TR/TE 2.7/1.1 ms, with 200 ms 
readout imaging duration. For heart rates above 
90 bpm a matrix of 192 x 130 was used. 



Images were reconstructed using a phase sensitive 
method which incorporates non-rigid correction of 
respiratory motion [22]. Co-registration between pre- 
and post-contrast acquired image series was performed 
in the generation of ECV maps [13]. Tl and associated 
SD-maps were generated using the PSIR MOCO method 
[22] and were processed in-line on the scanner, while 
Rl, ECV, and SD-maps of ECV with co-registration were 
processed off-line. 



In-vivo studies 

This study was approved by the local Institutional 
Review Boards of the National Heart, Lung, and Blood 
Institute and Suburban Hospital, and all subjects gave 
written informed consent to participate. Hematocrit 
used in calculation of ECV was measured from a venous 
blood sample drawn just prior to the CMR study. Tl- 
maps were typically acquired at least 15 min following 
administration of Gd contrast (0.15 mmol/kg) (Gadavist, 
Bayer Healthcare). In-vivo data was acquired to (1) meas- 
ure typical myocardial SNR, (2) to demonstrate the value 
of SD maps in assessing the statistical significance of sub- 
tle changes in Tl, and (3) to illustrate the application to 
measurement of ECV mapping. 
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Figure 3 Measured (dots) vs calculated (circles) standard deviation of Tl vs Tl (left column) for various image SNRs for 2 sampling 
strategies, 5(3)3 and 4(1)3(1)2, from Monte-Carlo simulation with 65 7 536 trials, and corresponding standard deviations for Rl = 1/T1 
(right column). 



Kellman et al. Journal of Cardiovascular Magnetic Resonance 2013, 15:56 
http://www.jcmr-online.eom/content/15/1/56 



Page 6 of 12 



0.07 
0.06 



■ 1 1 1 


Robust Fit Off / No outliers 

Robust Fit On / No outliers 

Robust Fit Off / single outlier 

Robust Fit On / single outlier 







-300 -200 -100 0 ICC 200 303 

T1 error (ms) 

Figure 4 Histogram of T1 errors for case of T1 = 1000 ms, SNR = 30 illustrating the difference between iteratively re-weighted robust 
fitting in the presence of outliers. 



SNR measurements 

The SNR in the myocardium was measured in-vivo in 
n = 20 subjects prior to contrast in order to derive typ- 
ical values of SNR used in the analysis. Images were 
reconstructed with scale in SNR units [29] and values 
were measured in manually traced septal and lateral wall 
ROIs. Reconstruction in SNR units [29] facilitates SNR 
measurement but requires off-line reconstruction from 
the raw acquired data. SNR maps were calculated from 
the signal intensities of the longest inversion time image 
which has achieved steady state. The MOLLI imaging 
protocol used for SNR measurement was modified to 
acquire 12 heartbeats to ensure complete recovery 
was achieved. The sampling for this protocol was 
12(1)1, i.e, images acquired for 12 heartbeats on the 
first inversion, followed by 1 recovery beat, and 
acquisition of 1 additional heartbeat on the second 
inversion. 



Results 

SNR measurements 

Myocardial SNR with the MOLLI protocol was mea- 
sured in 20 subjects using a voxel size of 1.4 x 1.9 x 6 
mm 3 . The SNR in the septum was 43 ±11 (m ± SD) 
and in the lateral wall was 22.8 ± 4.3. The SNR 
of the LV blood pool averaged over the full cavity 
was 42.7 ± 7.9. 



Monte Carlo simulation 

The measured and calculated estimates for Tl standard 
deviation as well as Rl = 1/T1 are graphed (Figure 3) for 
a range of Tl and SNR = 20, 30, and 40 which were in 
the range of SNRs measured using the specified imaging 
protocol. The 5(3)3 sampling strategy is simulated for 
Tl in the range 400-1600 ms, and the 4(1)3(1)2 sam- 
pling strategy is simulated for Tl in the range 300-500 
ms. The estimated SD is within 1 ms of the calculated 
SD in all cases, i.e., the estimate is unbiased. The vari- 
ability or standard deviation of the SD estimate was 
found to be related to the number of inversion time 
measurements used in the fit. For the 5(3)3 protocol 
with 8 measurements the standard deviation of the SD 
estimate was 1/3 of the mean SD value (mean over 
ROI), whereas for the 4(1)3(1)2 protocol with 9 mea- 
surements, the standard deviation of the SD estimate 
was 1/3.3 of the mean SD value. This variability in the 
distribution of errors was constant over the range of 
Tl s and SNRs simulated. This variability corresponds to 
a single pixel and is reduced by sqrt(N) when averaging 
over a ROI with N independent samples which was veri- 
fied in the Monte-Carlo simulation by averaging over 
blocks of variable length N. 

To illustrate the robustness to outliers (Figure 4), the 
histogram of Tl measurements for a single protocol 
(Tl = 1100 ms, SNR = 20, 5(3)3 sampling) was per- 
formed with a single randomly located outlier at 6 
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Figure 5 Tl-map (LEFT) for CuS04 doped agar gel phantoms for varying Tl and 12, corresponding SD map calculated from 200 trials 
(CENTER), and mean of estimated SD maps for 200 trials using proposed method (RIGHT). 



Kellman et al. Journal of Cardiovascular Magnetic Resonance 2013, 15:56 
http://www.jcmr-online.eom/content/15/1/56 



Page 7 of 1 2 



30 














25 




# / 




a 

(0 

H 15 

-0 

3 
en 

E 10 
% 

5 




* 




0 

c 




■ ■iii 






5 10 15 20 25 30 

calculated T1 SD (ms) 


Figure 6 Comparison of estimated (proposed method) versus 
calculated (200 repeated trials) SD for phantom data shown in 
Figure 5. 



standard deviations (n = 262,144 trials). Without outliers 
both the non-robust (plotted black) and robust (plotted 
blue) fitting have similar performance (SD = 30.3 and 
31.3 ms, respectively). With a single random outlier the 
SD becomes 71.6 ms (plotted red) vs 56 ms for the itera- 
tively re- weighted robust fitting (plotted magenta). 

Phantom measurements 

The proposed method of estimating SD was validated in 
agar gel phantoms using the method of repeated trials 
(Figure 5). The estimated values for SD correlated ex- 
tremely well (r = 0.997) with the calculated values across 
a range of SD s for tubes with varying Tl, T2 and SNR 
(Figure 6). 



In-vivo examples 

An example Tl-map for a normal subject (Figure 7) 
illustrates the variation in Tl estimation error with sur- 
face coil intensity variation. The anteroseptal ROI has 
SNR = 32.1, Tl = 1012 ms, SD = 25.0 ms, whereas the 
lateral wall ROI has lower SNR = 20.9, Tl = 1026 ms, 
SD = 41.8 ms (increased). The SNR in the LV blood pool 
ROI had Tl = 1580 ms, SNR = 29, and SD = 38 ms. 

Example Tl-maps illustrate detection of subtle, focal 
elevation in Tl in subjects with myocarditis (Figure 8) 
and HCM (Figure 9). The subject with myocarditis 
(Figure 8) has a focal elevation of the sub-epicardial 
region of the lateral wall (1098 ms) corresponding to 
103 ms elevation with respect to the septum (995 ms). 
The SD of in the lateral wall region is estimated to be 
43 ms, i.e., the focal elevation is approximately 2.4 SD 
on a pixel-wise basis indicating statistical significance 
(P < 0.01) for this individual which has even greater sig- 
nificance on a ROI basis (see Discussion). 

The subject with HCM (Figure 9) has a focal elevation 
of the septal region (1170 ms) corresponding to 84 ms 
elevation with respect to the septum (1086 ms). The SD 
in the septal region is estimated to be 36 ms, i.e., the 
focal elevation is approximately 2.33 SD on a pixel-wise 
basis indicating statistical significance (P ~ 0.01) assum- 
ing normal errors. The extension of error maps from Tl 
to Rl (Figure 10) and ECV (Figure 11) is illustrated for 
this subject with HCM. ECV and SD maps (Figure 11) 
are calculated from the AR1 maps after co-registration 
of pre- and post-contrast acquired maps [13] and scaling 
for HCT and AR1 of the blood pool. The ECV value in a 
septal ROI was 34.75% with SD estimated to be 1.0% 
and the lateral wall had ECV = 26.2 with SD = 1.2%. 
These SD values represent the pixel-wise variation. 

Non-rigid motion correction used in these studies is 
designed to correct respiratory motion and is not as 
effective in correcting cardiac motion due to significant 
variation in the RR interval. In these cases, motion 
related changes in the signal will lead to increased error 
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particularly around structures. Uncompensated motion 
related errors are readily apparent in the SD map 
(Figure 12) and may be used to assess the quality of 
the Tl-map. For the case with more severe cardiac 
motion (Figure 12 right) the subject had a regular pat- 
tern of premature ventricular contractions (PVCs). 



Protocol Comparisons 

The Tl-mapping protocols 5(3)3 and 4(1)3(1)2 used in 
this study are compared in Figure 13 with the widely 
used original MOLLI protocol 3(3)3(3)5 and the short- 
ened MOLLI protocol 5(1)1(1)1 [12] over a range of Tl 
at an SNR = 25. For the range of Tl values between 300 



T1 (ms) SD (ms) 




Figure 9 Example Tl and SD maps for subject with HCM acquired pre-contrast (TOP) and post-contrast (BOTTOM). Pre-contrast T1 maps 
exhibit focal T1 abnormalities in the septal region corresponding to T1 elevation of 84 ms relative to the lateral wall representing an elevation of 
2.3 SD on a pixel-wise basis (septal SD = 36 ms). 
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and 600 ms corresponding to subjects following adminis- 
tration of contrast, the 3(3)3(3)5, 4(1)3(1)2, and 5(1)1(1) 
are reasonably similar in performance. For longer 
values of Tl where it has been shown that the original 
MOLLI becomes heart rate dependent [11] the 5(3)3 



and shortened MOLLI have been proposed to mitigate 
this. The shortened MOLLI reduces the dependence on 
heart rate by discarding samples depending on the mea- 
sured RR interval. For the case of HR = 60 bpm shown 
in Figure 13, the precision of the Sh-MOLLI method is 



ECV (%) SD{%) 




Figure 11 ECV and SD maps for subject with HCM acquired corresponding to the Rl-maps shown in Figure 10, calibrated using LV 
blood and HCT values. 
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Figure 12 Example SD maps for illustrating cases with increased error due to uncorrected motion of varying severity. Uncorrected 
motion tends to cause local errors in the SD map which are readily recognized due to appearance of anatomical structure. The SD maps serve as 
an indication of T1-map quality. 



equivalent to the 5-0 sampling strategy for Tl > 1000 ms, 
and is equivalent to 5-1 between 400 and 1000 ms. 



precision of Tl and ECV estimates, and permits optimi- 
zation and comparison of protocols. 



Discussion 

Usefulness of the standard deviation measure 

The clinician is frequently faced with the question of 
whether a regional variation or a global shift is real or is 
a random fluctuation due to noise. The proposed gener- 
ation of an error map in meaningful units serves as an 
additional metric to aid in measuring statistical significance, 
as shown in the examples of Figures 8 and 9. As always, it 
is important to use judgement since there are other impor- 
tant sources of bias errors and artifacts. For example, 
regional variation in off-resonance can cause variation in 
apparent Tl which are treated here as bias errors since they 
do not increase the random fluctuation. This issue is 
important and must be considered in the assessment. The 
random component limiting precision is only a single 
component to consider, but serves as a performance bound. 
The proposed formulation may be used for predicting the 




Original 

MOLLI 

Protocol 

3-3-5 



Figure 13 Comparison of the precision of various sampling 
strategies employed in widely used Tl -mapping protocols. 



Accuracy vs Precision 

The issue of accuracy versus precision is an important 
subject. Although precision is essential to detecting sub- 
tle variations, the ability to detect global shifts in Tl is 
also limited by biases unless the normal reference values 
are reproducible, known and established. However, one 
must exercise caution since there are variable factors 
that may influence these "biases" such as off-resonance 
and shim, or protocol parameters such as matrix size. 
To some extent, the ECV may be less sensitive to accur- 
acy in cases where the percentage error is the same for 
both myocardium and blood for both pre- and post- 
contrast acquisitions; in this case errors could cancel. In 
order to realize the full potential of Tl and ECV as 
biomarkers, bias errors should be minimized and the 
magnitude of any bias variation should be known. 

SD estimation 

Estimates of SD were found to be unbiased when using 
the MAD approach for estimating the standard deviation 
of fit residuals. The use of iteratively re-weighted fitting 
to be robust in the presence of outliers comes at a small 
penalty of a few percent in statistical efficiency in the 
case for which the data is normally distributed (see dif- 
ference between blue and black plots in Figure 4 corre- 
sponding to 3% difference in SD). This sacrifice is 
deemed to be a small premium to pay for the gain in 
protection against outliers. 

The largest source of errors in pixel wise Tl -mapping 
is due to uncompensated heart motion due to variation 
in the RR interval, and this may be evident in the SD 
map. This increase in estimated SD due to motion error 
is apparent from structured appearance of the noise 
(Figure 12) and in this instance one must exercise 
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caution in interpreting the Tl-map values in these pixel 
regions which may be in error. 

Although the formulation has been developed and val- 
idated for PSIR reconstructed images, this formulation 
may be extended to saturation recovery (SR) methods in 
cases where the SNR is reasonably high such that the 
noise distribution is sufficiently normal. It is not readily 
translated to magnitude IR Tl -mapping due to the 
approximately Rician noise distribution in magnitude 
signal, particularly near the signal nulls. Furthermore, a 
frequently used implementation of the MOLLI method 
originally described [9,10] that uses the magnitude re- 
constructed images is based on a multi-fitting approach 
that perform multiple fits based on successfully in- 
crementing the unknown time of zero-crossing and 
choosing the value with the best fit. This procedure is es- 
sentially estimating a 4-th parameter (the zero-crossing) 
which degrades the precision for specific values of Tl 
which have zero- crossings near the inversion times sam- 
pled [22]. 

The error in measuring Tl was found to be well ap- 
proximated as a normal distribution, and therefore the 
analytic calculation agreed well with the Monte-Carlo 
results. 

Pixel-wise Tl and SD mapping 

Tl -measurement on a pixel-wise basis provides the abil- 
ity to detect regional heterogeneity and small focal eleva- 
tions that might be difficult to identify using ROI based 
fitting. It also provides a context for identifying boundar- 
ies between blood and myocardium to mitigate partial 
volume contamination of ROI regions. Although from a 
noise standpoint the precision of fitting to ROIs is better 
than the pixel wise fit, the precision of the pixel-wise fits 
averaged over an ROI are similarly improved. In the in- 
vivo examples shown in Figures 8 and 9, the ROI sizes 
in the region of elevated Tl were 48 and 150 pixels, 
respectively. The SD of Tl on a ROI basis is improved 
by the square root of the number of independent pixels 
averaged. The number of independent pixels used to cal- 
culate the improvement in SD due to the averaging is 
somewhat less due to partial Fourier and parallel im- 
aging factors and is estimated to be 40%, corresponding 
to 19 and 60 pixel averages, respectively. Thus the pixel- 
wise SD in Figure 8 example of 43 ms is approximately 
10 ms on a ROI basis, further improving statistical 
significance. 

SD maps for Tl and ECV may be used to provide confi- 
dence for gauging the statistical significance of abnormally 
elevated Tl and ECV on a pixel- wise basis, as well as indi- 
cating overall quality of the Tl-map. Although the SD 
maps for current protocols are noisy in appearance, meas- 
urement of the SD in typical sized ROIs are dramatically 
improved in addition to the improvement in the actual Tl 



values. The ROI size for the septal ROI in the Figure 9 ex- 
ample was 150 pixels. Assuming that the actual number of 
independent samples was approximately 60 due to inter- 
polated reconstruction, the SD estimate of 36 ms had a 
variability of 36/3/sqrt(60) = 1.5 ms. 

Conclusions 

Pixel -wise estimates of Tl mapping errors have been 
formulated and validated, and the formulation has been 
extended to ECV. The ability to quantify the measurement 
error has potential to determine the statistical significance 
of subtle abnormalities that arise due to diffuse disease 
processes involving fibrosis and/or edema and is useful 
both as a confidence metric for overall quality, and in 
optimization and comparison of imaging protocols. The 
formulation provided may be extended to other paramet- 
ric mapping measurements such as T2 or T2*. 
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